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Abstract 

The single diffractive production of the standard model Higgs boson is computed using the 
diffractive factorization formalism, taking into account a parametrization for the Pomeron structure 
function provided by the HI Collaboration. We compute the cross sections at next-to-leading order 
accuracy for the gluon fusion process, which includes QCD and electroweak corrections. The gap 
survival probability (< \ S\ 2 >) is also introduced to account for the rescattering corrections due to 
spectator particles present in the interaction, and to this end we compare two different models for 
the survival factor. The diffractive ratios are predicted for proton-proton collisions at the Tevatron 
and the LHC for the Higgs boson mass of Mh = 120 GeV. Therefore, our results provide updated 
estimations for the diffractive ratios of the single diffractive production of the Higgs boson in the 
Tevatron and LHC kinematical regimes. 

PACS numbers: 13.60.Hb, 12.38.Bx, 12.40.Nn, 13.85.Ni, 14.40.Gx 

Keywords: Higgs boson production, single diffractive, gap survival probability, next-to-leading order accu- 
racy, electroweak corrections 



1 



I. INTRODUCTION 



Many hard diffractive and high-energy processes are under intense study in the last 
two decades. One of the main subjects in high-energy physics is the improvement of our 
knowledge about quantum chromodynamics (QCD). Additionally, the Higgs mechanism is 
one of most important subjects to be investigated at the LHC, being a cornerstone in the 
electroweak sector of the standard model (SM). The Higgs boson is expected to be produced 
by the gluon fusion process (gg — > H), making the data analysis of this process an important 
topic in the project for the LHC experiments, like ATLAS and CMS. 

Recent analyses presented an updated estimation of the mass range where it is expected to 
observe the Higgs boson, which, combining the data coming from CDF and DO experiments 
at the Tevatron, have excluded the range 158 < Mh < 175 GeV with 95% of confidence 
level jl| . Furthermore, very recent simulations of the ATLAS experimental group have shown 
that a wider mass range can be excluded with the future LHC data. For instance, with an 
integrated luminosity of 2 fb _1 and a 8 TeV beam energy the range 114 < Mh < 500 GeV 
can be excluded with 95% confidence level 

The diffractive processes are well described by the Regge theory, where it is considered 
that a family of resonances is exchanged by the colliding protons |3j . The leading pole that 
accounts for this interaction will drive the high-energy behavior of the total cross section, 
being particularly labeled Pomeron, that has the vacuum quantum numbers [4|. However, 
the nature of the Pomeron is not completely known, as well as its reaction mechanisms, 

n 

but it is a successful formalism to describe hard diffraction data [5J. Moreover, based on 
the parton model, it was proposed that the Pomeron could have a partonic content, i.e., 
quarks and gluons as its constituents, by the Ingelman-Schlein (IS) formalism [6J. Then, 
systematical observations of diffractive deep inelastic scattering at HERA have increased 
the knowledge about the Pomeron, providing a diffractive distribution of singlet quarks and 
gluons into the Pomeron as well as the diffractive structure functions 3]. 

In this work we are interested in the single diffractive (SD) processes, characterized by the 
emission of a Pomeron from one of the colliding hadrons that scatters off the other hadron. 
The cross sections for the SD process are computed at next-to-leading order (NLO) accuracy 
with QCD and electroweak (EW) corrections, and we use the gap survival probability (GSP) 
from two different models that accounts for the survival factor for the diffractive Higgs 



boson production. The cross sections and the diffractive ratios are estimated for the process 
p + p(p) — >■ p + H + [LRG] + for the kinematical regime of the Tevatron (a/s = 1.96 
TeV) and for those expected to be reached in the LHC (i/i = 7, 8 and 14 TeV). In this 
approach the hard processes will occur by the interaction of the content of one hadron and 
the content of the Pomeron. In other words, the diffractive cross section is the convolution of 
the diffractive parton distribution functions (DPDF) and the corresponding partonic cross 
section, in a similar way as the inclusive case. In addition, diffractive events with a large 
momentum transfer are also characterized by the absence of hadronic energy in a certain 
angular regions of the final state, the so-called rapidity gaps. So, the SD processes will 
present in the final state a large rapidity gap between one proton and the Higgs boson as 
its main signature. 

For the Tevatron kinematical regime, it is known that the data are not correctly pre- 
dicted with the use of the IS formalism [8|, however there are important contributions from 
unitarity effects to the single-Pomeron exchange cross section that can be considered. These 
absorptive (unitarity) corrections take into account the fraction of large rapidity gap pro- 
cesses, except elastic scattering, being quite important for the reliability of the predictions 
for hard diffractive processes. The multi-Pomeron contributions depend on the particular 
hard process, and one is able to compute the GSP [9] for a specific production process, 
which accounts for the fraction of events where the rapidity gaps will be present in the 
final state after the rescattering events. In this way, the application of a survival factor in 



the diffractive cross section can correctly describe the high-energy data, 
predictions for W ± ) Z°, heavy quarks, T and J/i/j were presented in Refs. 



or instance, some 
for the LHC 



energies, and it was possible to see that this approach describes very well the Tevatron data. 

This paper is organized as follows: in Sec. [Til we present the main equations for the 
inclusive production of the Higgs boson at NLO accuracy. Next, in Sec. IHH we rewrite the 
parton luminosity in order to introduce the Pomeron exchange from the colliding proton, 
taking into account the gg — >• H production. Further, in Sec. IIVI we present the models for 
the GSP applied in this work, showing the probabilities for each energy regime. Then, in Sec. 
fVT. we present the estimations for the inclusive and diffractive cross sections as a function 
of the Higgs boson mass for different collider energies, and also the rapidity distributions of 
the Higgs boson. Finally, in Sec. IVI| we summarize our conclusions. 



II. INCLUSIVE PRODUCTION 



Let us present the main formulas for the inclusive cross sections for the production of 
Higgs boson in proton-proton collisions. The production process considered in this work 
is the gluon fusion pp — )■ gg — > H, since it is the leading production mechanism of the 
Higgs bosons in the high-energy regime The gluon coupling to the SM Higgs boson 



is mediated by a triangular loop of quarks, with the leading contribution of the quark top. 



The production cross section at lowest order is given by 
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with the Drell-Yan variable defined as th = Mfj/s, where s is the invariant pp collider 
energy squared. The gluon-gluon luminosity has the form 

d£99 f 1 dx 



-g(x : M 2 )g(r/x,M 2 ), 



(2) 



dr j T x 

with g(x, M 2 ) being the gluon distribution function into the proton, where we apply the 
MSTW2008 parametrization at NLO accuracy for such distribution [3], with M as the 
factorization scale. In Eq.(CE]), the function o"o reads 
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where A q (t q ) = 2[t q + (tq - l)/(r Q )]/r^, and t q = Mjj/Am 2 . In this work it is considered 
only the lea ding contribution of the top quark (m q = m t = 172.5 GeV), called heavy-quark 



limit in Ref. 



12J, and then we are taking the approximation tq < 1, which means the use of 



f(r Q ) 



arcsm 



The NLO QCD corrections for the 



usion process gg — >• H correspond to the processes 



gg — > H(g), gq — > Hq and qq — > Hg jl2,Q, introducing virtual and real corrections to the 
scattering amplitude. The production cross section for the Higgs boson at NLO accuracy in 



pp collisions is written as 
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with the renormalization scale in the strong coupling constant a s and the factorization scale 
in the parton densities to be fixed properly. Particularly, in Eq.© the strong coupling 



constant is applied at lowest order accuracy; however, for the NLO contributions the a s is 



applied at NLO accuracy through the exact numerical solution 
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din fi 2 An s lQn 2 
where (3 = (11N C - 2N F )/3 and ft = (102iV c - 38iV F )/3, with N c = 3. The A scale is fixed 
by the threshold of the quark masses during the /i 2 evolution, and fixing the value of Np 
properly. 

The coefficient C(tq) denotes the contributions from two-loop virtual corrections, reg- 
ularized by the infrared singular part of the cross section for real gluon emission, and is 
expressed by [12] 



,'llJV c -2JVp. \ . /x 
C(t q ) =n 2 + c(t q ) + ( c - *- ) log 
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(6) 



where n 2 refers to the infrared part, and c(tq) is a finite function, which, solved analytically, 
results in c(r Q ) = 11/2 for r Q = M 2 H /Am 2 q < 1 [16]. 

The AcTjj are the hard contributions from gluon radiation in the gg scattering and the qq 
annihilation, and they depend on the renormalization scale /i and the factorization scale M 
in the parton densities. These contributions can be expressed by 12 1 
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where f = tjj/t, and P gg (f) and P gq (f) are the standard Altarelli-Parisi functions [l 
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The F + denotes the usual + distribution, such that F(f) + = F(f) — 5(1 — f) L dr'F(r'). 
As we are considering the heavy-quark limit, the c?y functions can be solved analytically, 



resulting in a simpler set of expressions 12] 
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Finally, also included are the electroweak two-loop corrections 18], which enhance the 
total cross section by 5% in comparison to the NNLO QCD cross section. In this way, the 
total cross section is computed with the addition of the EW corrections by 

0"NLO = CQCD+EW = C"QCd(1 + few)- (10) 

The total cross sections for the inclusive process are shown by the solid curves in the Figs. 
[TJH] for different collider energies. The gray bands around these curves express the variation 
of the renormalization and the factorization scales in the range 0.5M# < (/i = M) < A.OMh- 
Looking particularly to the results for the LHC, our results reproduce the values obtained 



in Ref . 



12 



18], although it is not the case for Ref. [19 I 



III. DIFFRACTIVE PRODUCTION 



For the diffractive process, the calculations are based on the IS formalism for diffractive 
hard scattering p. In this case, the Pomeron structure is taken into account by its quark and 
gluon content through the parametrization of the DPDF. The SD cross section is assumed 
to factorize into the Pomeron-hadron cross section and the Pomeron flux factor. In other 
words, it consists of three steps: first, a hard Pomeron is emitted from one of the protons 
in a small momentum transfer \t\, being this hadron detected in the final state; then, the 
second hadron scatters off the emitted Pomeron; during the Pomeron-hadron interaction, 
partons from the Pomeron interact with partons of the hadron, producing the Higgs boson. 
Accordingly, we will take into account absorptive effects, multiplying the diffractive cross 
section by a specific survival factor for each collider energy. The luminosity for the SD 

18] and the curve presented in the Fig. 6 in Ref. 



1 Comparing the results obtained in the Ref. 



12 



19f for the 



total cross section in inclusive process, one can see that there is a disagreement between the results at 
y/ii=U TeV, since the NNLO cross section for M H = 200 GeV in Ref.[l9| 
predicted in Ref.flQ at NLO. 



is clearly smaller than that 
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The Pomeron structure function Fi/ W / P is expressed by 

Fi/w/p = fw/p(xjp)fi/ w ( — , M 2 J , (12) 

V^ip / 

with fw/p{zw) being the Pomeron flux, and fi/ w ((3, /i 2 ) the parton distribution into the 

Pomeron, where % stands for g, q, and q. 

In the estimates for the cross sections in Eq. ffTTj) . we consider a standard Pomeron 

flux from Regge phenomenology, which is constrained from the experimental analysis of 

the diffractive structure function. In this case, we apply the flux obtained with the HI 
II 

parametrization [7|. The Pomeron structure function has been modeled in terms of a light 
flavor singlet distribution S(x), i.e., the u, d and s quarks with their respective antiquarks. 
Also, it has a gluon distribution g(z), with z being the longitudinal momentum fraction of 
the parton in the hard subprocess. The gluon density is a constant at the starting evolution 
scale Qq = 2.5 GeV 2 . In our numerical calculations, we apply the cut x < Xjp < 0.05 in 
agreement with the HI parametrization. The Pomeron trajectory is assumed to be linear, 
a w (t) = aip(0) + a' w t, with «j and their uncertainties obtained from fits to HI forward 
proton spectrometer (FPS) data 20|. We choose x w f t mm fw/ P ^t = 1 at x w = 0.003, where 
\t m in\ ~ m p x ip/(l ~ x ip) i s ^ e minimum kinematically accessible value of \t\, m p is the pro- 
ton mass, and |t cu t| = 1.0 GeV 2 is the limit of the measurement. The HI parametrization 
provides two different inputs for the fit of the partonic structure functions. As our curves 
show very close results using both fits, we chose the fit A to perform our predictions in this 
work. 



IV. GAP SURVIVAL PROBABILITY 



In the diffractive cross sections [Eq.l fTT]) ]. we are further including the GSP < |>S| 2 >, being 



described in terms of absorptive corrections js]. It can be estimated using the equation 

2 f\A(s,b)\ 2 e-"^^b 

<|S|> " S\M»M>>?b ' (13) 



where A is the amplitude of the particular process of interest at the center-of-mass energy 
squared s described in the impact parameter space b. The quantity Q is the opacity (or 
optical density) of the interaction of the incoming hadrons. This suppression factor of a 
hard process accompanied by a rapidity gap does not depend only on the probability of 
the initial state survival, but it is also sensitive to the spatial distribution of partons inside 
the incoming hadrons, i.e., on the dynamics of the whole diffractive part of the scattering 
matrix. 

There are distinct approaches in the literature to compute the value of the < |5*| 2 >, pre- 
dicting different probabilities for the diffractive Higgs boson production. Applying a survival 
factor to diffractive processes brings an uncertainty to the predictions for the production 
cross sections jsj, since there is no accurate prediction for the GSP, resulting in an imprecise 
predictions. Hence, we compare two different models for the GSP, being the most applied in 
other works, in order to investigate the available calculations of the survival factor to drive 
our predictions, and certainly the ones that will be studied to describe the future data. The 



first one is that of Refs. [2l|] (labeled KKMR), which considers a two-channel eikonal model 
that embodies pion-loop insertions in the Pomeron trajectory, diffractive dissociation and 
rescattering effects. Then, the survival probability is computed for single, central and double 
diffractive processes at different collider energies, assuming that the spatial distribution in 
impact parameter space is driven by the slope B of the Pomeron-proton vertex. We will 
consider the value < |5'| 2 >kkmr = 6% (10%) for the SD process in the LHC (Tevatron). 



The second estimation for the survival factor is the model presented in Ref. 22] (labeled 
GLM), with a calculation for an eikonal single-channel approach. We take the case where 
the soft input is obtained directly from the measured values of atot, ®ei and hard radius 
Rh- The F1C approach was chosen to perform our predictions, resulting in a probability of 
< I 5 I 2 >glm = 8 -!% (12.6%) for the LHC (Tevatron) energy. We quote Ref. Q for a de- 
tailed comparison between this approach and the Kaidalov-Khoze-Martin-Ryskin (KKMR) 
one, including further discussions on model dependence of inputs and consideration of mul- 
tichannel calculations. 

Unfortunately, these models only account for the GSP in the kinematical regime of the 
Tevatron or the LHC energies, i.e., for yfs = 1.8 TeV and 14 TeV. In order to make precise 
estimations with reliable values for the GSP, we chose to adopt a similar way to estimate 

8 



the survival factor, in % for the desired energy, following the approach of Ref. 23| 

< i s ' 2> < % )-F^b- (14 > 

with the parameters a = 46.52 (30.77) and b = -3.80 (-4.41) for the GLM (KKMR) model. 
Then, the Table UJ summarizes all the survival factors applied in the predictions for the SD 
process. 

These particular models were chosen in order to indicate the uncertainty (model depen- 
dence) of the soft interaction effects. It is worth to mention that some implementations of 
GLM model include the results of a two- or three-channel calculation for < \S\ 2 >, which are 



considerably smaller than the one-channel approach 221 ] . 



V. RESULTS AND COMMENTS 

In this work we are mainly interested in the analysis of the cross sections for the SD Higgs 
boson production for different collider energies, bringing a rapidity gap in the final state as its 
main signature. Furthermore, the diffractive production can be an alternative way to detect 
;he Higgs boson in hadron colliders, since it provides a higher signal-to-background ratio 



24j . Moreover, the SD production cross sections are presented in Figs. [THU being the results 
presented with no survival factor by the dashed curves, and including the GLM (dot-dashed) 
and KKMR models (double-dot-dashed). Additionally, some values of the production cross 
section are presented in Table [Til for selected Higgs masses, showing specifically the values 
for the cross section with the adopted survival probabilities. As one can see, the production 
cross section in the kinematical regime of the Tevatron is very small, as expected. However, 
for higher energies the cross section reaches values of the order of 100 fb, showing that it 
may be possible to detect the Higgs boson in the LHC through the SD process. Besides, 
there are some detectors to be set up at the LHC experiments, and they will make it possible 
to tag the outgoing proton [25] or to detect forward showers |26j. Then, the SD events can 
be an effective way to look for the Higgs boson at the LHC. 

Still, to have a clear analysis of these results and to estimate the fraction of diffractive 
events, we compute the diffractive ratio for the Higgs boson mass of Mh = 120 GeV, being 
presented in Table IIII1 As expected, the diffractive ratios are small and growing with the 
collider energy. Specifically, the GLM model shows a ratio nearly constant from the Tevatron 



energy until 7 TeV, and then growing for higher energies. However, the KKMR model shows 
a different behavior, presenting a decrease in the same region, and growing slowly at higher 
energies. This effect is observed due to our assumption for the survival factor for the collider 
energy of 7 and 8 TeV. However, a proper calculation of the survival factor for these energies 
may show a higher GSP than those presented in the Table [H which will increase the ratio, 
reaching a similar behavior as the result for the Tevatron. 

In fact, these values show that the SD events will have a very small rate in the LHC 
kinematical regime for a luminosity of a few fb _1 . It encourages the implementation of 
specific detectors in order to detect the rapidity gaps or the forward protons, since the 
Higgs boson discovery from its decay products is going to be more difficult in the inclusive 
production due to the huge background signal [27]. 

Finally, in Figs. 5 and 6 we present the rapidity distribution of the Higgs boson for 
different collider energies. For higher energies, the distributions are clearly central, which 
shows that the contributions from the parton distribution function and the DPDF have 
larger values in central rapidity. Particularly, as the momentum fractions of the parton 
into the hadron A increases, the one of the parton into the hadron B (or the Pomeron) 
decreases uniformly, achieving a higher combined contribution for yn = 0. However, this 
same behavior does not occur in the Tevatron kinematical regime, showing two distinct 
peaks in the distribution where the parton distribution function and DPDF have its higher 
combined contribution. In the results for mid energies (7 and 8 TeV), one can see that the 
distributions are still central, however showing very small peaks around \yn\ = 2. 



VI. CONCLUSIONS 

In summary, we have evaluated estimations for the SD Higgs boson production in pp 
collisions at the Tevatron and the LHC, considering the IS formalism with the introduction 
of rescattering corrections. We are using the Regge factorization to calculate the SD cross 
sections at NLO accuracy (QCD+EW). In particular, we take a parametrization from HI 
Collaboration for the Pomeron structure function, extracted from their measurements of 
F^ 3 \ with the results directly dependent on the quark singlet and gluon content of the 
Pomeron. For the available fits in this parametrization, we chose the fit A to perform our 
predictions. For instance, the cross sections are of about cr SD = 0.4 (0.1) pb for M H = 120 

10 



GeV for yfs = 14 TeV (7 TeV). These cross sections are higher than that obtained 
77 production mechanism, predicting a production cross section of 0.12 - 0.18 fb 



xom the 



28|], and 



even for the exclusive Higgs boson production 221]. Moreover, comparing our estimations 
with no survival factor, we pre dict a cross section for = 14 TeV higher than the previous 
results for the SD process [121]. In addition, the two models considered for the GSP have 
computed a survival factor that has a variation of 25%. This difference is significative in 
order to perform reliable predictions for the Higgs boson production; however, it is expected 
that they are going to be tuned with the future data from the LHC experiments. In any case, 
our predictions with different survival factors may give a good estimation for the production 
cross section for the presented collider energies. 

The SD production of the Higgs boson may not provide significative advantage in com- 
parison to the inclusive production, since the background can not be suppressed using the 
same statements as in the double Pomeron exchange case. Nevertheless, the hadronic activ- 
ity in the final state will be reduced in the SD events, increasing the possibility of observing 
the Higgs boson. Furthermore, the rapidity gaps may be observed in the LHC with the use 
of specific detectors, and then it can bring new data to be compared to the SD estimations. 
These results are the first NLO predictions for the single diffractive Higgs boson production, 
applying the GSP in the diffractive factorization. Besides, we have feasible values for the 
diffractive cross sections, and diffractive ratios as well, but the difference in the predictions 
is a bit high, which reveals that a study of the GSP for the multiple-Pomeron interactions 
in SD events is highly necessary. Therefore, we have presented updated estimations for the 
diffractive Higgs boson production, allowing the possibility to compare them to the future 
LHC data. 
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<|5| 2 > (%) 


GLM KKMR 


1.96 


12.3 9.7 
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9.2 6.9 
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8.9 6.7 


14. 


8.1 6.0 



TABLE I. Estimations for the survival factor for different collider energies obtained with Eq, (|14l) . 



Mass (GeV) 


(TeV) 






1.96 7. 8. 


14. 


120 


5.36(4.23) 88.59(66.44) 119.70(90.11) 


346.43(256.62) 


140 


2.57(2.02) 58.69(44.02) 81.43(61.30) 


248.75(184.26) 


160 


1.24(0.98) 39.56(29.67) 56.07(42.21) 


183.06(135.60) 


180 


0.60(0.47) 27.60(20.70) 40.23(30.28) 


134.46(99.60) 


200 


0.31(0.24) 19.96(14.97) 29.10(21.90) 


104.65(77.52) 



TABLE II. Estimations for the production cross section (fb) at NLO accuracy for selected Higgs 
masses. The values are shown for the GLM (KKMR) models for the GSP and Mh = /jl = M. 



Ratio (%) (TeV) 

1.96 7. 8. 14. 
R S d 6.23 8.31 9.10 11.21 
Rglm 0.76 0.76 0.81 0.90 
Rkkmr 0.60 0.57 0.61 0.67 

TABLE III. Estimations for the diffractive ratios in different collider energies for a Higgs boson of 
Mjj = 120 GeV. The ratio is predicted for single diffractive events (-Rsd) an d for both models of 
the GSP (-Rglm and -Rkkmr) with Mh = (jl = M. 



15 




I I I I I I I I I 

100 120 140 160 180 

Higgs mass (GeV) 



FIG. 1. Production cross sections for inclusive and single diffractive processes for the Higgs boson 
in the LHC. The gray bands show the variation of each cross section for the energy scales 0.5Mjj < 
(// = M) < A.OMh- The lower curves show the predictions for the single diffractive production 
using two different models for the GSP. 
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FIG. 2. The same as Figfflfor y/s = 8 TeV. 
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FIG. 3. The same as Figfflfor y/s = 7 TeV. 
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FIG. 4. The same as Fig(T]for the Tevatron energy. 
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FIG. 5. Rapidity distribution (pb) in single diffractive process of the Higgs boson for the current 
energy of the LHC (7 TeV) and for its future kinematical regimes. 
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FIG. 6. Rapidity distribution (fb) in single diffractive process of the Higgs boson for the energy of 
the Tevatron (1.96 TeV). 
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